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Abstract 



We measure the mass of the r lepton to be 1775.1 db 1.6(stat.) ± 1.0(sys.) MeV using 
r pairs from Z° decays. To test CPT invariance we compare the masses of the positively 
and negatively charged r leptons. The relative mass difference is found to be smaller than 
3.0 x 10~ 3 at the 90 % confidence level. 
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1 Introduction 



We have measured the mass of the r lepton, m T , using data taken by the OPAL detector during 
LEP running at the Z° resonance. The first test of CPT invariance using r leptons is performed 
by measuring the masses of the positively and negatively charged r leptons separately. To 
determine the r mass we use two pseudomass techniques, previously established by ARGUS 
and CLEO 0, which rely on the reconstruction of the mass, energy, and direction of the 
hadronic system in hadronic r decays. 

In the OPAL detector charged particle tracks are reconstructed by a central detector consist- 
ing of a silicon microvertex detector, a vertex chamber, a large jet chamber, and z-chambers. 
Photons coming from the decay of neutral pions are measured using a hermetic lead glass 
calorimeter located outside of the solenoid and an iron-streamer tube sandwich calorimeter is 
used to measure hadronic showers. This hadronic calorimeter is used in conjunction with muon 
chambers mounted around it to separate muons from hadrons in the event selection and decay 
mode identification. The OPAL detector is described in detail elsewhere [J3J. 



2 The Pseudomass Method 



In a hadronic r decay, the mass of the r lepton is related to the 4-momentum of the resulting 
hadronic system^ by the formula 



m 2 ,, 



(Pr-Ph) 2 = rr%-2E T E h + 2\p T \\p h \ cos + ml, \p T \ = ^/E 2 T -m 2 T (1) 



where p T (ph), Pt (Ph), and E T (E h ) are the 4-momentum, 3- momentum, and energy of the 
r (hadronic system). The neutrino mass, m u , is set to zero and the r energy is taken to be 
the beam energy, thus permitting the mass of the r lepton to be reconstructed if the angle ip 
between the direction of the r and the hadronic system were known. Due to the high boost 
at LEP energies, tp is limited to a few degrees, with the exact limit on ip depending on m^. 
A pseudomass, m* n (one value per cone|), is defined by taking ip to be zero. It gives the true 
mass when ip = and is smaller than the true mass in any other case. The distribution of m* cn 
for r — > 3 7r ± v T decays, presented in fig. [I] is a broad distribution with a sharp cutoff at the r 
mass. Further, the details of the shape of the distribution are determined by the dynamics of 
the decay with the general tendency that the more massive the hadronic system, the closer the 
pseudomass is to the cutoff. However, the position of the cutoff only depends on the mass of 
the t lepton. The r mass is extracted from this position. From fig. [I] a small tail of pseudomass 
above the cutoff is evident. It is caused by three effects: background, resolution, and initial 
state radiation, the latter weakens the assumption E T = E^ CSim . More details can be found 
in M. 



^^Here 'hadronic system' is defined as the sum of all the hadrons produced in the decay. The sum of their 
momenta gives the momentum of the hadronic system. 

2 Cones with a half-opening angle of 35° arc defined around the leading particles in the event. In most cases 
a cone contains the decay products of exactly one t. 
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Figure 1: The distribution of the pseudomass m* n from the r — > 3 7r ± z/ r class. The points 
with error bars are data and the open histogram is the Monte Carlo prediction including the 
background from misidentified r decays (shaded area) normalized to the data. The solid line 
shows the parameterization in the fit window (1.6 to 2.0 GeV) after the final unbinned maximum 
likelihood fit to the data. The events have been binned for the presentation only. 

In r pair events where both r leptons decay hadronically, a second pseudomass, rra* v , can 
be defined, which is derived from the event as a whole, not just from one of the two cones. 
In such events the acollinearity a, i.e. the angle between the hadronic systems from the two 
decays, carries additional information on the r mass (see fig. |2[). This acollinearity restricts the 
decay angles of the r leptons through the relation + ip_ + a > it, where ip + and ip- are the 
decay angles of the positive and negative r in the laboratory frame, respectively. If the sum of 
these three angles were known, the mass of the r lepton could be reconstructed on an event by 
event basis. Since this information is unavailable the pseudomass, m* v , is calculated instead, 
assuming that the sum of the three angles is 180 degrees and that the masses of the positive 
and negative r are equal. The m* v distribution shows similar features to m* n . An example is 
shown in fig. [| Again, the cutoff at high values, from which the r mass is extracted, is clearly 
visible. More details can be found in ||. 



h 




Figure 2: Illustration of the relevant angles: ip + and tp_ are the opening angles between the 
r lepton and its hadronic decay products, a is the angle between the hadrons of the two r 
leptons. In the absence of radiation, the r leptons are produced collinear. 
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Figure 3: The distribution of the pseudomass m* v from the class of r — > n v T decays recoiling 
against r — > tt^tt ^. The points with error bars are data and the open histogram is the 
Monte Carlo prediction including the background from misidentified r decays (shaded area) 
normalized to the data. The solid line shows the parameterization in the fit window (1.5 to 
2.0 GeV) after the final unbinned maximum likelihood fit to the data. The events have been 
binned for the presentation only. 



3 Selection and Classification 



The selection of r pair events is performed as described in [|J using the data recorded during 
the years 1990 to 1995 at center-of-mass energies at and around the Z° mass. A total of 159 373 
t pairs are identified within a geometrical acceptance limited to polar angles^ of | cos^l < 0.95. 
The overall efficiency including the geometrical acceptance is 89 %. The purity of the sample 
is 98 %. The background consists mainly of electron and muon pairs. 

Each selected event is split into two narrow cones of particles, one for each r. All 1-prong and 
3-prong conesQ with a minimum momentum of 5 % of the beam momentum are then subjected 
to a likelihood classification in order to identify the decay mode of the two r leptons using the 
same procedure as in [0] . Table [I] summarizes the decay modes considered. 

For the calculation of the pseudomasses m* n and m* v , all tracks are assumed to origi- 
nate from charged pions. Neutral pions are reconstructed from showers in the electromagnetic 
calorimeter with the algorithm described in [||]. The number of neutral pions reconstructed in a 
cone does not necessarily match the number of 7r° expected from the decay mode identified by 
the likelihood selection. For example, there might only be a single 7r° reconstructed in a cone 
identified as a r — > 7r ± 2 7r° v t decay. In such a case the missing neutral pions are ignored in the 
calculation of the 4-momentum of the hadronic system. If there are too many ir° reconstructed, 
we exclude those with the lowest energies. 

3 We use a righthanded coordinate system with the z-axis pointing in the direction of the electron beam and 
the x-axis pointing towards the center of LEP. 

4 A 1- or 3-prong cone is a cone containing 1 or 3 charged tracks. 
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1-prong modes 


3-prong modes 




eff. % 


pur. % 


decays 




eff. % 


pur. % 


decays 


r — > TT^ 1 v T 

T — ► 7T^~ 71^ V T 

T — > 7T^~ 2 7T^ V T 


76.2 % 
60.2 % 
45.1 % 


64.5 % 
77.1 % 
40.9 % 


41148 
58518 
26093 


r — > 3 n^- v T 

T — > 3 7r^~ 7T^ z/ r 


71.5 % 

61.6 % 


83.8 % 
64.2 % 


24832 
11636 



Table 1: Efficiency, purity, and number of decays selected by the likelihood identification in 
the decay modes explored in this analysis. In addition to these modes, the likelihood identifies 
leptonic decays and various modes with one charged hadron plus additional tracks from photon 
conversions or Dalitz decays. All numbers refer to the full range of pseudomasses. The purity 
in the region of the cutoff is substantially better. 
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Table 2: The assignment of the various decay modes from the likelihood identification into 
classes analyzed by either m* n (cn) or m* v (ev). The row labeled 'no ID' combines the cones 
where no identification was possible with those identified as leptonic decays and all other 
channels not used in the analysis, r — > 371^ v T includes r — > ^ K° v T . For r — > tt v t versus 
t — ► 3 7T^~ u r see text. 



The pseudomass used to analyze an event depends on its identified decay modes. This 
choice of method can either be m* n applied to one or both cones of the event independently 
or m* v applied to the whole event. Table illustrates the association of the decay modes of 
the r to the two methods. The association has been optimized to give the smallest possible 
statistical error on m T , a process which results in decays with high hadronic masses typically 
being assigned to the m* n method, leaving the lighter hadronic masses for m* v . Every decay or 
event is only analyzed by one method to avoid statistical correlations. 

A slightly more complicated procedure is used for events with a r — > 3 7r ± v t decay recoiling 
against a r -> 7r ± v t . The 3-prong cone is initially subjected to the m* n method and the 1- 
prong cone neglected (see Tab. 0). If m* n turns out to be below the fit window described in 
the following section, the event is analyzed through m* v instead. This procedure could also 
be applied to other combinations of decay modes, but the gain in overall resolution is only 
marginal. 
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class 


method 


entries 


result in 


MeV 








1467 


1777.6 ± 


7.2 


7r±2 7T° 




™cn 


1311 


1801.7 ± 


12.6 


3 7T ± 




<n 


2680 


1776.1 ± 


1.9 


3 TT+TT 




™c*n 


3467 


1777.6 ± 


6.2 


71" ^ <->• 


7T ± 


m e * v 


483 


1790.2 ± 


9.4 


71" ^ <-» 




m e * v 


1346 


1759.9 ± 


6.4 


TT^ <-> 


7T ± 2 7r° 


m e * v 


506 


1779.9 ± 


16.4 


7T ± 7r° <■ 


-> 7T ± 7r° 


m* ev 


1191 


1770.6 ± 


8.7 


Tt^TT <- 


-> 3 7T ± 7r 


Kv 


565 


1793.3 ± 


13.5 


7T =t: <-> 


3 7T ± 


m* cv 


328 


1757.9 ± 


10.9 



Table 3: The 10 classes of the analysis. The method applied, number of entries in the fit 
window, and a result (statistical error only) from a fit to each class alone are specified. The 
'7r =l= <-> 3 7r ± ' class only contains those events not used in the '3 7r ± ' class (see text). 

In the following steps of the analysis the distributions of m* n are handled in exactly the 
same way as those of m* v . We define the following classes (see Tab. 0): 

• Four classes corresponding to the decay modes r — > n° u T , r — > ^ 2 7r° u t , t — > 3 v T , 
and t — > 3 7r ± tt° v t which are analyzed with the m* n method; 

• Five classes for the different combinations of decay modes analyzed with the m* v method; 

• One class for the r — > v T versus r — > 3 7r ± v t decays described in the previous para- 
graph. 

Each class is further subdivided into subclasses according to the expected resolution. This 
avoids having information from well reconstructed entries degraded by the overlay of those 
entries with poorer resolution. There is a total of 25 subclasses, between one and five per 
class. The criteria used for the classification of the expected resolution is the quality of the 
reconstruction of the hadronic system, features like the presence of z-chamber or silicon hits 
on the tracks, the number and momentum of the reconstructed neutral pions, the probability 
of a vertex fit or the polar angle of the cone. The best resolutionF] in the pseudomass m* n of 
17 MeV is achieved in the r — > 3 7r =l= v T class, and this class also has the largest weight in the 
final result. For m* v , the best resolution is 21 MeV in the r — > v T versus r — > v T class. 
Table |3| gives some more information on the classes. Figures |T] and ^ correspond to the subclass 
with the best resolution. 

5 The resolution is denned as the r.m.s of the difference between true and reconstructed pseudomass from 
simulated events that fall into the fit window. 
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4 Extraction of the r Mass 



The value of the r mass is extracted from the pseudomass distribution by means of unbinned 
maximum likelihood fits. The distributions are parameterized by a step function multiplied 
by a polynomial of up to third order. For each subclass one of the following step functions is 
chosen, depending on which gives the best description of the data: 



r i 1 ' 



/(*) 



2 V VTT 



l + e x 

1 /TT 



(2) 



, — arctan [x / 

TT V2 



The variable x is related to the pseudomass by x = hi (a, + 5 + m*) where a, and foj are fit 
parameters and m* is the pseudomass, the measured quantity. The index % labels the sub- 
classes. The parameter 5 is a common shift of all distributions. The order of the polynomial is 
increased from zero until a reasonably good parameterization of the distribution is achieved. Q 
Only events in the vicinity of the cutoff are used for the fits. The upper end of the fitting 
window is 2 GeV for all subclasses. The values of the lower bound range from 1.4 to 1.6 GeV, 
depending on the subclass. For each subclass the lower bound is decreased from 1.6 GeV until 
no further significant reduction of the expected statistical error is observed. After the initial 
parameterization of the m* distributions, the fit windows are reduced by 50 MeV on each side 
to ensure a stable parameterization in all areas that might be reached in the course of the 
following fits. The type of step function, the order of the polynomial, and the size of the fitting 
window are determined for each subclass separately from a Monte Carlo simulation ||. 

The functional dependence of the pseudomass distribution on m T is very simple in the 
vicinity of the cutoff. In this region the shape of the distribution does not depend on m T . On 
a variation of m T the whole distribution is shifted along the pseudomass axis by an amount 
identical to the shift in m T . We will use this relation to extract the result, although it is not 
strictly correct away from the cutoff. We have checked that the fit windows are small enough 
so as not to introduce a systematic bias from this procedure. 

The r mass is extracted in three steps. In the first step, the distributions from simulated 
events are fitted for each subclass separately in order to fix the parameters of the fit function, 
i.e. the position and width of the step function (aj, and the coefficients of the polynomials. 
In the second step the same simulated distributions are simultaneously fitted using the slightly 
reduced fit window. A common shift <5mc is the only free parameter in this fit with all other 
parameters fixed to the values obtained in the first step. In the third step, the same common 
fit is applied to the data, but with 5a a ta as the free parameter. The result is calculated from 
m T = m™ c + cWa - S M c with c = 1777 MeV the r mass used in the simulation. The result 
is 1775.1 MeV with a statistical error of 0.4 MeV from the Monte Carlo and 1.6 MeV from the 
data. Results were also determined for each class separately as a cross check, the results are 
shown in table 0. 



5 We histogram the data and apply a \ 2 test to decide on the quality of the description. 
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5 Systematic Uncertainties 



The most significant contributions to the overall systematic uncertainty are described in this 
section. 

• The dominant systematic uncertainty comes from the calibration of the tracking chambers 
and the electromagnetic calorimeter. We use muons and electrons of 45 GeV from Z° 
decays and Bhabha scattering to check these calibrations. We find uncertainties in the 
momentum scale of the tracking chambers relative to the Monte Carlo of 0.1 % and for 
the calorimeter of less than 0.3 %. Several different scenarios are employed to extrapolate 
these uncertainties to the momentum and energy ranges relevant to r decays. The most 
pessimistic scenario, a scaling of 1/pt by a factor, gives a systematic error from the 
tracking of 900 keV. The uncertainty in the calibration of the calorimeter is found to 
have a negligible impact (250 keV), since it only effects channels with neutral pions. 

• Uncertainties in modeling the resolution of the tracking chambers result in an error of 
40 keV on the measurement of m T . The resolution of the electromagnetic calorimeter has 
a negligible impact. 

• In order to estimate the systematic error contribution from the modeling of the dynamics 
of the r decays, we have varied the mass and width of the ai meson in the 3^ final state 
by ±100 MeV and changed the fraction of ujtx events in the 47r final state by ±20 %. The 
first two variations each change the result by 100 keV, while the third has a negligible 
impact. We also tested two different models for the r — ► 3 7r ± ^ T decay 0. Using these 
models the result changes by less than 100 keV. 

• In the derivation of the pseudomass formulae the assumption E T = E hca _ m is used, which 
is only true in the absence of initial state radiation. Low energy photons from initial state 
radiation have only a negligible effect on the measurement, while a variation of ±1 % of 
the rate of hard initial state radiation results in an uncertainty on the measurement of 
70 keV. 

• The uncertainty in the calibration of the beam energy of LEP has a negligible effect. An 
error of 2 MeV on the beam energy changes the result by less than ±40 keV. 

• The calculation of the pseudomasses assumes that the v T has zero mass. A value of 
18.2 MeV, the current limit on the r neutrino mass ||, increases the result by 110 keV. 
No systematic error is assigned. 

• Systematic uncertainties related to background from misidentified r decays are negligible, 
since their pseudomass distributions show no significant structure in the region of the 
cutoff. Background from non-r sources is also negligible. 

• To address possible biases from the parameterization of the pseudomass distributions, we 
increased the lower edge of the window for all channels, in steps of 25 MeV, to 1.6 GeV. 
Furthermore, we fitted the distributions with the type of step function giving the second 
best description of the data. Both studies showed consistent results. 

Adding in quadrature all systematic uncertainties including the statistical error from the Monte 
Carlo gives a total systematic uncertainty on the measurement of 1.0 MeV. 
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6 CPT Test 



To test CPT invariance we compare the masses of the positively and negatively charged r 
leptons. The pseudomass method m* n allows a separate measurement for the positive and the 
negative r lepton, whereas the m* v method implicitly assumes the two masses to be identical 
and therefore the m* n method is used for all decays. The measurement does not use any Monte 
Carlo simulation which reduces the systematic uncertainties. 

For the CPT analysis the m* n distributions are separated according to the net charges of 
the cones. To extract the result we use a procedure similar to the one used to measure m T . 
In the first step, the distributions of the positive or negative r leptons are fitted subclass by 
subclass to fix the parameterizations. Then all the distributions of the positive r leptons are 
fitted simultaneously with the parameters fixed to the values from the first fit, allowing only 
for a common shift in mass 5 + . Then negative r leptons are fitted the same way with S~ as the 
free parameter. If CPT invariance is conserved, the two shifts must be equal. We obtain the 
result 5 + - 5~ = 0.0 ± 3.0 MeV. 

In a r pair event produced in e + e~ collisions, a mass difference between the positive and 
negative r will also create a difference in energy between the two r. Although in principle this 
invalidates the assumption E T = E heam , this effect is numerically negligible. 

Most sources of systematic errors affect the result for the positive and negative lepton in 
the same way, so that their contributions cancel. The largest error in the mass difference 
comes from possible differences in the calibration between positively and negatively charged 
tracks which is limited to less than 0.2 % by studying muon pairs. This implies a systematic 
uncertainty on the mass difference of 1 MeV. All other systematic errors are negligible and the 
total systematic error is small compared to the statistical error. It is added in quadrature to 
the statistical error. 

The m* n method uses the momentum and mass of the hadronic system measured from the 
curvature of the tracks in the magnetic field of the detector. These tracks originate mainly from 
pions with a small fraction coming from kaons. The charge and the charge to mass ratio of 
these particles have to be known in order to convert the observed curvature into a momentum 
measurement. However, without assuming CPT invariance, it is no longer obvious that the 
charge and mass are the same for positively and negatively charged pionsQ. Experimentally, 
the relative pion mass difference has been measured to be (2 ± 5) x 10~ 4 which is below our 
sensitivity, but the pion charge difference has not been measured separately. However, the pion 
charges could be different due to a charge difference between the r + and the r _ or because 
of charge non-conservation. A charge difference between tc + and n~ as large as 5 % would 
invalidate our measurement. 

This is the first test of CPT invariance with r leptons. The result is 



where m(r) is the charge-independent result from the previous section. This result implies that 
the relative r mass difference is smaller than 3.0 x 10 -3 at the 90 % confidence level. 

common deviation of the charge of positive and negative pion from the charge of the electron or the r 
would not spoil the measurement. Only a difference between 7r + and tt~ is relevant. 




(3) 
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7 Summary 



We have measured the mass of the r lepton from the pseudomass distributions of the hadronic 
r decays recorded by the OPAL detector at LEP and obtained a result of 



1775.1 ± 1.6(stat.) ± 1.0(sys.) MeV 



(4) 



where the first error is statistical and the second systematic. The result is compared to other 
measurements in fig. |j. Using the pseudomass m* n we obtain an independent measurement 
of the mass of the positive and negative r. The two numbers are consistent and a limit of 
3.0 x 10~ 3 on the relative mass difference is be placed at the 90 % confidence level. This is the 
first test of CPT invariance in r physics. 
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Figure 4: Comparison of our result with the results from ARGUS [|TJ, CLEO 0, and BES |10 
The shaded band indicates the current world average 
result. 



which does not include the OPAL 
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